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Abstract. We provide a uniform decay estimate for the local energy of gen- 
eral solutions to the inhomogeneous wave equation on a Schwarzchild back- 
ground. Our estimate implies that such solutions have asymptotic behavior 
|0| = O^r~-^(|t — |r*|| 2^ as long as the source term is bounded in the norm 

{'^-^)~^ ■{l + t+\r*\)^'^L^{H^{r'^dr*duj)). In particular this gives scatter- 
ing at small amplitudes for non-linear scalar fields of the form Dgi/fi = X\4i\''4' 
for all 2 < p. 



1. Introduction 

In this paper our goal is to give a somewhat elementary discussion of the global 
decay properties of general solutions to the scalar wave equation on the exterior of 
a Schwarzchild black hole. That is, we consider the manifold with boundary: 

(1) M = Rx [2M,oo) X §2 ^ 
with Lorentzian metric: 

(2) ds - -(1 - — ) dt^ + (1 - —)-'dr' + r^dio^ , 

r r 

and we look at smooth functions (p which do not touch the boundary of |^ for each 
fixed value of the parameter t and which satisfy the inhomogeneous wave equation: 

□g0 ^ \7°'da(t) = G. 
The main question we would like to answer here is: How quickly does the local 
energy of the wave (j) dissipate over compact sets in the r coordinate, and how 
precisely does the dissipation depend on the source G? 

In the special case of Minkowski space, M — 0, a quite satisfactory answer to this 
question is known. Here one has the classical uniform local energy decay estimate 
of C. Morawetz: 

(3) 

' il + u^)\m^ + il + u^)\L(f>\^ + il + u^ + u^WcP\^+ ^ ' ' " I'^l' ) 

2 



3x{t} 

< 



([ Ul + \u\ + \u\)G{s)\\L2^d.)ds) + [ (l + r2)|Vt,,0|2 . 



The second author would like thank MSRI and Princeton University, where a portion of this 
research was conducted during the Fall of 2005. 
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Here one sets u = {t — r), u — (i + ?'), L ~ 29„, and L = 29„. For the original 
proof see the paper of Morawetz 10 , and for an alternative proof as well as many 
generalizations see the recent work 0. 

The beauty of the estimate Q is that it gives one a huge amount of information 
about the global dispersive properties of the function 4>. For one, it produces a 
pointwise in time decay of the Lf^^ norm as well as the local energy. What's more, 
this local decay is given in such a way that it is clear what is happening on the 
whole of each time slice t = const., even very far away from the origin r = 0. In 
fact, using Sobolev embeddings and only rotations, the Morawetz estimate is good 
enough to provide uniform decay at the rate of (1 + t)~^. However, perhaps the 
most important property of the estimate J^Jl is that it turns out to be incredibly 
useful when dealing with non-linear problems. This is because it places very simple 
conditions on the source term G, the kind which are relatively straight forward to 
recover in bootstrapping arguments given the form of the left hand side of Q . For 
instance, ((SJ makes dealing with the global existence problem for small amplitude 
non-linear scalar fields of the form: 

(4) 00 = X\cP\P<P , 

essentially trivial in the case where 2 < p. All one has to do is to combine (jSJ with 
appropriately localized Sobolev embeddings to yield the decay estimate: 

which is enough to integrate the nonlinearity on the right hand side of Q when it 
appears on the right hand side of Q. In fact, if one takes into account characteris- 
tic estimates of the form see again 0, then it is possible to push the exponent 
p to certain values p ^ 2. We will not discuss such refinements here. 

What we will show here is that for the more general case of M ^ 0, an es- 
timate which is essentially of the form (j3Jl holds in the case of general Schwartz 
(on each fixed time slice) functions 4>. The proof we give is a relatively straight 
forward integration by parts, similar in spirit to how one proves (O- In the final 
section of the paper, we indicate how our estimates can be used to give a short proof 
of global existence and decay for non-linear scalar fields of the form Q when 2 < p. 

To state our main theorem, we will use the following notation. We first reparametrize 
the radial variable in the usual way: 

(5) r* ^ r + 2M ln(r - 2M) , 

and then introduce the optical functions and null-generators for the coordinates 
{t,r*,Lo): 

(6a) u = (t + r*) , u = [t~r*) , 

(6b) L — dt + dr* , L = dt — dr' ■ 

We will prove that: 

Theorem 1.1 (Uniform Local Energy Decay for the Scalar Wave Equation on 
Schwarzchild Space). Let (i,r*,w) he the coordinates (as defined above) on the 
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manifold M. Q with metric ((5J) . Let (f) be a smooth function compactly supported 
on each hypersurface t — const, and set: 

Og^ ^ G . 

Then one has the following global estimate: 



(7) 



xS2x{t} 



{l + u^){L{rcp)f + {l + u^){L{rcp)f 



(1+^2+^2^(1 ) 

r 



1 M 



dr*cLj 



< 



f /■* 2Af \ ' 

^ II (1 + \u\ + \u\)(\ - _) r • (v/1 - A,p„G)(s) \\L-(ir'i^) dsj 



XS2x{0} 



2M . ( 1 



M 



dr*duj 



where the implicit constant depends only on the mass M . Here Asph is the Laplacean 
in the angular variable lo, and)^^ is the associated gradient} 



Remark 1.2. Let us first give a heuristic summary of the content of the estimate 
{7)) and how it contrasts to the situation of flat space M = 0. Roughly speaking, 
the Schwarzchild space can be split into three pieces where one sees quahtatively 
different behavior in solutions to the wave equation. 

The first region is very close to the boundary r — 2M. For the static space-time 
we are considering, this is quite easy to understand. Here wave propagation looks 
essentially trivial in that one has (j) ^ Fi{t + r* ,u!) for some smooth decaying func- 
tion Fi on the space R x S^. That is, wave propagation near the boundary r = 2M 
is essentially just transport in the variable u = t + r*. The caveat is that this 
variable is the only one in which it is possible to get decay for this region because 
Fi{0,uj) does not need to be small. 

We would like to call the readers attention to the fact that the precise decay of 
the function Fi in the first variable seems to be quite delicate issue, and we will only 
obtain \Fi{u,uj)\ < (l+|u|)-5. Of course, this is all one should expect given that 
the right hand side of is consistent with initial data decaying at this rate. Since 
this estimate does not ask for a lot of information, which is actually its strength in 
treating non-linear problems, it does not give a lot of information in return. For a 
much more precise asymptotic in the case of spherical symmetry, and for the more 
difficult case of dynamic space-times (for the spherically symmetric coupled scalar 
field), we refer the reader to the very deep recent work of Dafermos-Rodnianski (see 
0-0) 



^This should not be confused with y from hue which is the covariant derivative on spheres 
of radius r. Of course these two only differ by the factor of r~^. 
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The second region is in the "far exterior" t < r* where things look essentially 
flat. This is also fairly easy to understand. Here one expects that things look very 
similar to the well known case of Minkowski space. 

The third region is "the boarder" close to r = 3Af , which in Regge- Wheeler 
coordinates ^ we extend to the region \r* \ ^ it. This is by far the most difflcult 
region to understand, and where one looses regularity in the estimate (UJ. This loss 
of regularity is in sharp contrast to the estimate © in the case of Minkowski space, 
and is also something one sees only by looking at the non-spherically symmetric 
(functions) situation. What is happening here is that the geometry is pulling the 
radiation apart into the two regions just mentioned, and there is a danger that this 
"splitting" could allow some fairly large residual portion of the radiation to linger 
for a long amount of time in the transition region r* ~ 0. 

Now it turns out that this effect can only happen (and it does happen!) if the 
wave 4> has a very high angular momentum. In this case it can concentrate on a 
very small set for each fixed time in the uj variable, and it will essentially rotate 
around the sphere r = 3M for a long while before dispersing. This behavior can 
be understood by observing that null geodesies tangent to the surface r — 3M will 
remain tangent to this surface [Hj, and that high angular momentum solutions to 
the wave equation will closely follow these geodesies for a long period of time before 
dispersing. 

This slow dispersion can also be understood by conformally changing the met- 
ric l(2Jl by the factor (1 — ^^). Since the coefficient of dt^ is constant on this 
new manifold, the corresponding wave equation describes the time evolution of a 
wave on a three dimensional Riemannian manifold with metric given by the spatial 
portion of the confromal Lorentzian metric. A simple calculation shows that this 
Riemannian manifold has a totally geodesic sphere (and hence closed geodesies) 
at the value r = 3M. Now, the original wave equation is equivalent to the wave 
equation with respect to the conformal metric modulo a smooth potential. For 
very high frequencies this potential cannot compete with the principle part of the 
conformal wave operator, so it is not difficult to construct coherent state solutions 
which concentrate near the closed geodesies sitting at r = 3M for a long amount of 
time. Therefore, from this point of view, one should look at (T)) as a sort of "cheap" 
dispersive estimate, and it is well known that such estimates loose regularity when 
the underlying geodesic flow is not well behaved in the sense of spreading of the 
classical trajectories. 

We further remark here that the nature of the surface and geodesies at r = 3M 
can be a little confusing to discuss in the relativistic terminology. The null geodesies 
at r = 3M which orbit the black hole form a helix in four dimensional space-time 
with an axis in the time direction. Although their projections onto the three dimen- 
sional coordinates (r, uj) is closed, because the t coordinate is constantly increasing 
they are not closed null geodesies. Also, while it is true that the hypersurface 
r = 3M is foliated by null geodesies it is not itself a null hypersurface, which is one 
with a null normal direction (see j^). Note that a normal to r = 3M is dr* which 
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is space like. 

As far as our analysis is concerned, the presence of null geodesies at r = 3M 
manifests itself through trapping terms which are positive for r ~ 3M. For a 
wave equation with a potential Q, we refer to x ■ \7Q + 2Q as the trapping term. 
This expression appears as a contribution governing the growth of the conformal 
(Morawetz) energy. It can be seen as the main "error" which is generated by the 
divergence of the conformal energy density, and is given by the first two terms on 
the right hand side of equation (|50|l below. This identity was first derived using a 
some what different formalism in the dissertation of the first author (see 0). 



Remark 1.3. Our proof of the bound Q will be very general in the sense that we 
derive it from a fairly generic family of estimates that holds for 1-D wave equa- 
tions with "strongly repulsive" potentials. It is to be hoped that this procedure 
can be used to accommodate other situations, such as higher spin equations on 
Schwarzchild space and possibly other space-times where spherical harmonic de- 
compositions make sense. We will leave these discussions to further work. 

The approach we take here is based on the previous works which proved 

space-time Morawetz type estimates on Schwarzschild-like manifolds, and the the- 
sis ^ which proved a version of the conformal (Morawetz) energy estimate ^ 
with growing right hand side. In the estimate contained in 01, the trapping term 
(described previously) appears integrated in space-time against the quantity t(<^)^, 
where (j) is the scalar field. If the factor of t were not present and if the field cj) 
were restricted to a single spherical harmonic, then the Morawetz estimate from 
would be sufficient to control this space-time integral. However, due to the fact 
that the reduction to individual spherical harmonics leads to trapping terms which 
grow quadratically according to the angular frequency, both an additional angular 
derivative and the factor of t must be controlled. 

In this paper, we present a simple argument which allows one to absorb the 
trapping term with the factor of t, and hence prove In the dissertation ^ and 
the forthcoming work ^ , a more involved phase space analysis is used to reduce the 
loss of angular regularity in the space-time Morawetz estimate and in the analogue 
of 121 to only e powers of the operator 1 — Agph- We leave the combination of these 
two techniques to future work. 

Acknowledgments: First and foremost we would each like to thank our advisors 
Matei Machedon and Avi Soffer for introducing us to the problem of decay estimates 
for the wave equation on Schwarzchild space, and for generously sharing their ideas 
with us. Without them this paper would, quite simply, not exist. We would also 
like to thank Hans Lindblad for reading a preliminary version of this article and 
giving several useful comments, as well providing a good deal of support and overall 
encouragement. Special thanks also to Jason Metcalfe, Igor Rodnianski, and Chris 
Sogge for many helpful conversations along the way. 
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2. Preliminary Setup 

In this section we will set up some preliminary notation and ideas from one 
dimensional wave equations on Minkowski space. This material is for the most part 
entirely standard, and we make no claim of originality for the basic concepts. Now, 
it turns out that Theorem is actually a special case of a family of estimates 
which holds in the following general situation. We consider 1-D wave equations of 
the form: 

(8) - Q{x)^ = H , 

where □ — —df + d^ and Q{x) is some smooth real valued function which we assume 
is general for the time being. When the source term H vanishes the field ^ comes 
from a Lagrangian with energy momentum tensor: 

(9) T^pii^] = 5„va/3V - (a^MV- + Q(^) • W) • 

A quick calculation using the equation (jH)) shows that one has the divergence iden- 
tity: 

(10) d^T^fsm = -i9^(Q)-(^)2 + i7.a^V . 
Also, the trace identity: 

g"^T^M = -Q{x) ■ {4,f . 
follows immediately, where g = diag{—l, 1) is the 1-D Minkowski metric. 



The utility of the tensor © is that 
to the flow of various vector-fields X = 
momentum density associated to X: 

(11) (^)P. 



is keeps track of how the field (O reacts 
X°'da on R X R. In general, we form the 



— TapX^ , 



and we compute from (|1U|) the divergence: 

(12) a" = -\x{Q) ■ i^)' + \t^p^''^t:^P + H ■ X(V') , 

where ^-^^i: is the deformation tensor: 

= d^X(, + dpX^ . 

In the next section we will use this setup to prove the following general 1-D uniform 
local energy decay estimate: 

Theorem 2.1 (1-D Morawetz Estimate for Positive Strongly Repulsive Potentials). 

Let tjj be a function on (1 -I- 1) Minkowski space which is compactly supported for 
each fixed value of the time variable t. Suppose ip satisfies the equation © for some 
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smooth real valued function Q{x) which satisfies all of the following conditions: 

(13) Q , (Positivity) 

(14) -xd,^{Q) (Repulsive 1) 

(15) xd.^{Q){x) ^ 2Q{x) -Csgn{x)d,^{Q){x) , x i Bi , (Repulsive 2) 

(16) xd^{Q){x) + 2Q(x) s; C\x\^^Q{x) , x B2 , (Homogeneity) 

(17) {l + X^)x^ s$ -Cxd:^{Q){x) , xe2Bi, (Critical Point) 

(18) s; Q{x) ^ C(l + A^) , xe 2Bi , (Local Bounds) 

where C and A are fixed non-negative constants with C strictly positive, and the 
Bi are compact intervals containing the origin. Then the following uniform local 
energy decay estimate of Morawetz type holds: 

(19) / ((l + M')(ii/^)' + (l + w')(^V^)' + (l+2i' + u')Q(x)-(V')^) dx 

jRxit} 

+ I l|(l + N + |u|)(l + A)i7(s)|U2(^,).||(l + A)(|Vt,,V| + Q^-|^|)(s)|U2(^,) ds , 
Jo 

+ [ Ul + h\ + \u\)H{s)U2^,,yEHi:{s))ds. 

Here we have set: 







U = t + X , U = t ^ X , 

L = dt + dx , L = dt- dx , 

and: 
(20) 

M^is)) ^ I {{I + u^){Li,f + [l + u''){L4,f + {I + 1? + u^)Q{x) ■ {ijf) dx 

JVlx{s} 

The implicit constant in estimate (|19|) depends only on the constant C and the size 
of the two intervals Bi, and not on the value of t or X or any other property of the 
potential Q{x) then those listed above. 



2.1. The case of Schwarzchild Space. Before moving on to prove the estimate 
l)19|l. let us first briefly indicate how this can be used to prove the bound In 
the {t,r*,uj) coordinates one writes the wave operator \g\^^ dag°'^\g\'^ dp as: 

il~^r'{~d^^ + r-^dr.ir'dr^^)) + ^A,phcl^ ^ G. 

Here Agph is the Laplacean in the lo variable. Introducing now the quantities xp = rcj) 
and H — {1 — ^^)rG this last line becomes: 

(21) -d^i, + d^,.i;-r-'d%{r)xp+^- j^A^p^^j = H . 
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We now follow the usual procedure of projecting this equation onto individual 
spherical harmonics. Since all of our estimate are both and rotation invariant, 
there is absolutely no harm in doing this. We write: 

where the form an orthonormal basis for the space lS.sphy\ — — A^Yx- On each 
harmonic the equation (|21() becomes: 

(22) -a,VA.. + a2.^;,,,-r-ia2.(r)^A.. -^^r^^\^ - Hx,^ . 

Dropping the (A, i) indices, labeling r* — x, and using the notation: 

xafi _ IK) 

(23) Qx - r-ie(r)+ \^ \ 
equation H22(l becomes: 

-d^i; + d^^-Qxix)^ = H . 

We now wish to apply the estimate to each of these equations, after we apply a 
spatial translation by a quantity a;o(A) which will be determined in a moment. The 
resulting family of estimates can then be safely added to obtain the full estimate 
as long as one can produce a single point xo{oo) such that xq{X) —^ xo{oo), 
and a single set of objects {C, 81,82) such that the conditions H13() - H18|l hold for 
(C, Bi + xo{X),B2 + xo(A)). Luckily, for the family of potentials (|23|l where ^ A 
is any real number, this is easy to show. The reader should keep in mind that the 
reason this is possible is that the conditions (|13|I - H18|I are not really size conditions 
on the potential Q\, but are actually conditions on the sign of Qx and its first de- 
rivative. This type of condition is very stable under multiplication by large positive 
constants, so it is not hard to get uniform behavior for large A. We will state this 
result as follows: 



Lemma 2.2. Let Q\ be the family of potentials defined on line (|23|) above, and 
set Xo{oo) ~ 3A/. Then there exists a constant C, a pair of sets Bi,B2, and a 
family of points a;o(A) — > a;o(oo) such that the potentials Q\ {x + a;o(A)) satisfy the 
conditions H13(l - H18|l for the triple (C,Bi + a;o(A), +xo(A)). All of these objects 
are completely determined by the value of M. 



Proof. First, notice that condition l|13|l is immediate. Next, recall that in the 
current notation we have x — r{x) + 2Af ln(r(a;) — 2M). We now write: 

^ , ^ , 2Af, /2M A2\ 

(24) Qa(x) = (l-_)(^_ + _j . 

The proof centers around showing that Q\ has an isolated critical point. We com- 
pute the first derivative: 

2M 2M /2M A2\ 2Af 2 /6Af 2\^ 

(25) = -4(1 - — ) ■ [A'r2 - 3A/(A2 - l)r - SAf^] 
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The polynomial on the right hand side of this last expression has exactly one root 
for positive values of r. This is given by the quadratic formula: 

(26) ^(A) ^ 3Af (A^ - 1) + Af v/9(A2 - 1)^ + i2\^ 

We now show that this positive root is trapped inside the interval [^^^ 3M]. Since 
it is clear from (|26|l that asymptotically r(A) — > 3Af, it suffices to show that r(A) 
is an increasing function for ^ A. This follows at once from differentiating the 
polynomial on line (|25|l with respect to the parameter A which yields the identity: 

_ 6AfAr-2Ar2 
' ^ 2A2r - 3A/(A2 - 1) ' 

A simple calculation shows that this quantity is indeed positive whenever r e 
[^M,2>M]. Therefore, we shall pick our sequence of points a;o(A) according to the 
rule r(a;o(A)) = r(A). This immediately gives the condition (|14|l for the family of 
translated potentials Q\{x + xq{X)). Also, note that the pointwise bound H18|) is 
immediate for any compact interval Bi. 

We now show the critical point behavior H17() . This boils down to the fact that 
the polynomial on line (|25(l has a simple root at r(A). In fact, taking the second 
derivative of the potential Q\ with respect to x and evaluating at the point xq{X) 
we have that: 

2 

Notice that this quantity never vanishes, and is 0{X^) as A oo, so one has H17|) 
for any compact set Bi given a suitable constant C, independent of the value of A. 

It remains for us to show the "strongly repulsive" conditions (|15|) - (I16|I hold for 
a uniform constant C and pair of sets Bi. This follows from direct inspection of 
the formulas (|24|l and H25() . We consider the cases of a; ^ ±oo separately. In the 
case of a; — > oo we also have that r oo, and we have the two asymptotic formulas 
(with uniform constants in A depending only on the mass M): 

1}? 6A/ A^ 1 
{x - xq{X)) ■ Q'x{x) = + O(^) + O(^) + {something negative} , 

„^ , , 4A/ 2A2 ^ 

^Q\\x) = — 5- ^ ^ + {somcthmg negative) . 

Notice that the {something negative} terms on the right hand side of the first line 
above contain logarithmic corrections of the form ln(x)/a;'* and A^ ln(a;)/x^, which 
come from the second summand on the right hand side of jSJ. It is important 
these come with a good sign. Now, combining these last two lines we have that as 
a; ^ 00: 

{x-xoiX))-Q'^ix) + 2Qx{x) ^ -^ + 0(^) + 0(-l) . 

This is enough to imply H15|) - H16|l because as a:: ^ 00 we also have the following 
strict lower bounds: 

^A*^ „, ^A^ 1 

73' 



^ ^ -Q' , V ^ -Q^i^) 



X 
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Finally, we deal with the bounds H15(l " (|16|l as a; ^ —oo. This is even easier to 
treat. Notice simply that both Q and d^{Q) are O ((1 + ~ ^)), while the 

factor (x — xq{X)) goes to — oo. This means that the first term on the left hand side 
of both (|15|l ~ (|l(j() is a very large negative multiple of the second. Therefore, the 
bounds (|15|I - H16|) are trivial because the left hand side is asymptotically negative. 
This completes our demonstration of Lemma 12.21 □ 

To wrap things up for this section, let us just mention briefly how one can sum 
the estimate (|19|) over the angular frequency localized pieces rpx.i- The key thing 
here is that the estimate (|19f) has been set up in such a way that one can use the 
Cauchy-Schwartz inequality to deal with the terms on the right hand side of (|19|l . 
Specifically, summing this bound over (A, i) indices and using the fact that the {Y^} 
form an orthonormal system on the sphere §^ we have that: 

(27) sup E{Hs)) = sup / ({l+y?)iLir<j))f + {l+u^){Lir(t>)? 

+ {l + u +u )■{ + ^{r<l)) )) dr duj 

< E [(l + A^)I(0)(V'v) 

+ J II {i+h\+W\){i+x)HxA^) hndr-yW {l+X){\Vt,r'^X,.\+{^-^)H^+^)^\^xMs) hndr-) ds 

+ / \\{^ + \u\ + H)Hx4s)\\L2(^dr'}-RH^x4s)) ds , 
Jo J 

where we are defining: 

M^xAs)) = I ({^ + u^){L4>x.^? + {I + n'){Mx.f 

JRx{s} ^ 

Now, bringing the sum under the integral sign in the two terms on the right hand 
side of H27|) above and then using the L^-L°° Holder inequality yields the bound: 

(28) 

sup M^is)) < ^((1-A,p^)^(0))+ sup [Ei(bis))+Ei{l-A,ph)his))]~' 

ft 2M 

■ II (1 + N + l"l)(l )r- Wl-A,phG){s) h^dr-du.) ds , 

Jo ^ 

where the usual energy is given by: 

JKxS^xjs} V r r r J 
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The estimate (0 now follows from ()28|l and taking an angular (momentum) deriv- 
ative of the basic energy estimate (see the next section for a proof): 

sup Ei^is)) < Eim) + if \\{^-—)rG\\L2i^dr'd.)ds 



3. Proof of the Main Estimate 

We now turn to the proof of the estimate (|2.1() . This will be accomplished in 
a series of three steps, each of which represents a tightening of the usual energy 
estimate. These are: 

(1) Usual conservation of energy. 

(2) Weak local decay of energy. 

(3) Strong uniform local decay of energy. 

Steps (1) and (3) involve a direct use of the energy-momentum tensor identities 
recorded in the previous section applied to various vector-fields X which are asso- 
ciated with the various types of decay as just listed. To prove item (2) above we 
use a Soffer-Morawetz type multiplier similar to what was done in 

Step 1: Conservation of energy. This is well known. In the current setup it 
comes from setting X — T — dt m Hll|) . Because Q{x) is time independent and 
since T is Killing we end up with an essentially divergence free quantity: 

C)"(^)P„ = H-dt{iP) . 

Integrating this over a time slab we arrive at the energy estimate: 

(29) / {(LiP)^ + {L^f + Q{x) ■ i^bf) dx 

jRx{t} 

< / II H{s) \\LHd.y\\ dMs) Wmdx) ds+ f {{L^f + {D^f + Q{x) ■ (V)') dx , 
Ja jRx{o} 

where the implicit constant is fixed and does not depend on Q (it is easy to calculate 
but we will not bother). 



Step 2: Weak Local Decay of Energy. In this subsection, we prove that the 
local norm of the quantity Q ■ (ip)'^ decays sufficiently fast in an average sense. 
Our bound will be rather weak in that we allow the right hand side of the estimate 
to grow like Xt. However, this weak bound will be precisely what we need in the 
next subsection when we prove the strong uniform local decay of energy. What we 
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propose to show is the foUowing estimate for integers 1 ^ iV: 



(30) 



/ / {l + s)Q-{i,fdxds - j j (l + s)xi{^)sgn{x)d.AQ)-{i^? dxds 




< sup N-'^E{s) + N{1 + X^)E{0) 



+ N 



JQ JR 




Here the imphcit constant depends only on the constants C and the lengths of the 
interval Bi from lines (|13|I - H18|) . We are defining E_ as the Morawetz type energy 
from line (|20|) above. Finally, xi denotes a smooth bump adapted to the interval 
[— 1, 1], and E denotes the basic energy from line (|29|l above. 

In our proof of l|3U|) it will be convenient for us to make the assumption that the 
local bound (|17|l drastically improves if we restrict to very small sets containing 
a; = 0. In particular, we will assume that: 



for some sufficiently small parameter e which will be chosen in a moment. It is 
crucial for us to point out here that our choice of e will only be dictated by the 
constant C and the size of Bi, and will not depend on any other property of Q. 
Also, it is immediate that the assumption 1)3 l|l in fact involves no loss of generality. 
This is because the equation © rescales as follows: 



i;{t,x) ^p{€^H,e-^x) , Q{x) e-^Q{e-^x) , H{x) e-^H{e~h,e-^x) . 



Notice that the conditions (|13|l - (|18|l adapt to the rescaled situation in obvious ways. 
In particular one has %M\ on the set Bi = e ■ Bi. For the rest of this subsection we 
will work in the rescaled situation where we assume all of H13|l -H18 |l as well as (IHlll . 
Of course once one has in this rescaled situation, one can recover the same 
bound for the original potential Q by scaling back. This will create constants which 
depend on e, but we choose this parameter only to overcome two things. The first 
is the possibly large constant C on the right hand side of H17|l (which is actually 
only a problem when A is small). The second is the fact that the original Bi may 
be small, so the constant cg^ on line lp?T|l where our improved bound holds is also 
small. Of course both C and cgj are fixed no matter how much we rescale, so these 
can be made up for by taking e sufficiently small. The main thing to keep in mind 
here is that our rescaling will never create constants in our estimates which depend 
in other ways on the shape of Q, other than the original assumptions we have made 



To prove Ij^UII . we use the following growth multiplier of Soffer-Morawetz type: 



(31) 



^■x^ ^ -xd^{Q){x) , 



x\ < CB^e , 



(nnj-iiHi). 



A(s, x)tl) 



(l + s)xi( 



10a; 
1 + s 
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where Lp is defined as follows: 
(32) ^{x) 



1 



dy , 



/o (1 + 12/1)'= 

where 1 < fc is a fixed constant. In practice the smaller the value of k the more 
favorable the estimates, so the reader may assume that k = 2. However, we will do 
all of our calculations in the general case so the overall structure is more apparent. 
The estimate H30II will follow from the usual procedure of directly calculating the 
integral: 



/ = — / / H ■ A{s , x)!^ dx ds , 
Jo Jr 

(dtijj - dlijj) ■ A{s, x)iIj dx ds + 



Qtp ■ A{s, x)ili dx ds 



(33) 



h 



and then using a Poincare type lemma near the critical point of Q{x) to get rid of 
the factor —xd^iQ) in favor of Q. We now compute the terms 1^ separately and in 
order. The first term Ii is the pure error. We first integrate by parts with respect 
to dt which yields the identity: 



(34) h = 



dtip ■ A{s){dtip) dxds 



dtij ■ A{s){ip) dx ds 



x{t} 



dtip-A{t){^) dx- 



x{0} 



dttp ■ A{0){ij) dx 



Here the operator A{s) is given by 

A{s)lp : 



, 10a; . 10a; , , 10a; 

xi(t-^) - T—^ii-r—) 



■ 1 + s 1 + s 
Also, one has the adjoint formula: 



1 + s' 



{l + \x\r 



lOr 

A*is)^ = -Ail;-20x[{-—)ip-^ . 

1 + s 

Therefore, a bound for the absolute value of the right hand side of (|34|l above is: 



(35) < n xi{^ 

Jo Jr 1 + s 



(dt^^f + (a^v-)' 



+ sup 



lOx 



(i + s)xi(t— )• (dt^py + {d^^^py + 

1 + s 



(l + |a;|)2'= 
1 



dx ds 



dx 



Ofis^tjRx{s} + * \ (l + |a;|)2* 

where xi is another [—1, 1] adapted smooth bump. To deal with the terms involving 
the inverse |a;| weight we use the Poincare type estimate: 



(36) 



1 



ii^y dx < ii^y{o) + 



.0 (i + N)' 

This follows at once from evaluation of the integral: 
1 



{dxi>y dx 



(i + |a;o|) 



((^)2(xo) + (V')2(-xo)) - 2(V)2(0) 



sgn{x)d^ 



(1 + W 



dx 
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and using the Cauchy-Schwartz inequality. Using now H36|l and the condition H18|l 
it is easy to bound: 

by using a partition of unity on Bi and M \ . 
Our next step is to use the bound: 

1 + s ^ ' 

where the right hand side is the Morawetz type energy from hnc (I20|) above. This 
and the bound (|35|l aUows us to estimate: 

< / (1 + s)~2 -^(s) + sup {^ + s)-^-E{s) 

Jn{1 + \) JV(l+A)sJs^t 
rN{l+\) 

E{s) ds + sup (1 + s) • E{s) , 

0^s<JV(l+A) 

< sup (A^(l + A))"^^(s) + sup 7V(1 + A)£;(s) . 

Using now the energy inequahty H29|l to deal with the second term on the right 
hand side of this last line we arrive at the bound: 

(38) < (l + A)-i(R.H.S.)(|2nj ■ 

In a moment we will need to multiply all of our estimates through by the factor 
(1 + A), so (|38|) is of the correct form. 

Before moving on to the second integral on line H33|l above, we mention briefly 
how to take care of the first integral on the right hand side immediately above that 
line. Applying the Cauchy-Schwartz inequality we have the bound: 



(39) / / \H{s)\-\A{s)^p\ dxds 

II (1 + s)H{s) |U.(,,) • II (|9,VI + XI (1 + Ixir^im^) hHd.) ds . 



< 

Using now a Poincare type estimate of the form (|37f) to deal with the last term on 
the right hand side of H39(l easily yields: 

(L.H.s.)13ni) < (1 + A)-^ (R.H.S.)!!! , 

which is sufficient for our purposes. 

Finally, we deal with the integral I2 on line 1)33(1 . After several integration by 
parts (this is an essentially well known calculation) we arrive at the identity: 



^2 = E ^: 



4 

3 ' 
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where the integrals Ki are: 



Ki = 

^3 = 



10x'i( 



10a; 



1 + s 

1 Or 



/-, X / IO2; , 



dx 



Bounding the first three terms above is essentially the same as what we have just 
done for the term Ii above. One simply uses Cauchy-Schwartz, the Poincare esti- 
mate H36(l ■ and the definitions of the two energies E and E_ to prove that: 

(40) \Ki\ + \K2\ + \K3\ < sup {N{l + X))-^E{s) + sup N{1 + \)E{s) . 

Therefore, the heart of the matter now is to obtain a positive lower bound for the 
quantity K4 in such a way that we can estimate the left hand side of H30|l . 

Before continuing with the proof, let us make a further simplification. Without 
loss of generality we may assume that the cutoff function Xi is the square of yet 
another smooth cutoff function, say xi- This allows us to replace "0 by xiV' in 
K4 above modulo a term involving [c^^jXi] = 0{j^) which is also cutoff where 
a;| < lO^^t. It is clear that this will again be of the form (1 + A)^^(R.H.S.)(|2ni), so 
we can just tack this error on to (|40|) above. 

Thus, what we will need to show is that there exists a sufficiently small constant 
c such that the following reverse bound holds for compactly supported functions tJj: 



(41) 



1 



2^'(9,^)2 - ^d^{Q) ■ (^)2 ^ ^"'(^) 



dx 



> c 



dx 



Once this is established, the bound (|30|1 will follow from combining the bounds 
l)38|l and H40|l with H41|) and the following estimate which also holds for smooth 
compactly supported functions ip: 



(42) 



Q ■ (i/jf dx 



Bi 



sgu{x)dx{Q) ■ {'4>f dx 



< (1 + A) 



dx 



We first prove H41|) . The overall strategy for this is very simple. The main thing 
we will establish is that the form of the weight function (|32|l reduces everything to 
having a "good" bound for the function tfj ai x = in terms of the hand side of 
l|lT|l . This latter task is relatively easy to accomplish because the assumption lp?T|l 
essentially means that —xdx{Q) ^ e~^i5o, where i5o is the unit mass at the origin. 
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This means that the potential term on the right hand side of (|41ll wiU give us a 
bound on ^p{0) with an 0(e2 ) constant. The details of this procedure are as follows. 
We first compute: 



= / 9, [if" dx , 

(43) = / if'" dx + 2 I if" ^pd^^ dx 



It will now be useful to have the identities: 

/,/ ^ -k ■ sgn(x) „ oi X : + 



Therefore, the right hand side of H43|l and a Cauchy-Schwartz gives us the bound: 



(ipy dx 



^ 2 



(l + |a;|)'=+2 



We may now assume without loss of generality that in this last bound we have 
C ^ B, otherwise there is nothing to prove on line H41|) . Therefore, dividing 
through by and squaring this last line we arrive at the bound: 



(44) 



ip'" i^f dx 



S - C s$ 4— ^A + 4 ( — I At-Ci 
k+1 \k+l 



The dangerous term is now the second one on the right hand side above. This needs 
to be controlled in terms of a sufficiently small constant. In fact, we will show that 
it is 0{e2) times the (R.H.S. 1 1)4111 . which implies that it may be safely absorbed 
into half of the remaining portion of A and a small amount of the potential term 
on fR. 11.8.1 141)1 . The bound which allows us to do this is the following: 



WW < e ( / T^^ZKk dx + I ~pd^{Q) ■ {iPf dx 



1 

(1+1^1) 

From the assumption ()31|l . this last estimate follows from: 

(45) (V)'(O) < e([ 2\ ..t. id.i'f dx + e-^ f x\{e-' x) ■ {i,f dx 

where x is some smooth 0(1) bump function whose support depends on the size 
of the set Bi from line 1)171) . The estimate 1)45)) is essentially scale invariant, so it 
suffices to show that: 

(46) (^)'(O) < / ^^—^id,^fdx+ ( x\{x) ■ {^f dx . 

This, in turn, follows from cutting things off and using the usual Sobolev embedding 
once we have the bound: 

(47) WxHWl^ < I j—^Y^{d..^fdx+ [ x\{x)-W'dx, 
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for some slightly smaller cutoff function x. This last bound can be proved in two 
steps. We first show the estimate: 

(48) m^ti^Wl. < [ —^-^{d:.^fdx+ [ x\{x) ■ {xPf dx , 

for some intermediate cutoff x- This bound follows at once from evaluating the 
integral: 

= / sgn(a::) 9^ x^ X' (V") 



dx , 

and using the Cauchy-Schwartz inequality to bound the error terms by (R.H.S.)(@SIl- 
Having now established (|48|l wc can prove (|47ll by applying the same procedure to 
the integral: 



= / d^[xx- (^)'J dx . 
Jr 

This completes our proof of (|45|l . and hence our demonstration of the main com- 
mutator estimate (|41f) . 

Having now dealt with the bound (|41|l , the only thing left for us to do in our proof 
of H30|l is to show the bound (|42|l . Notice that the bound for the second term on 
the left hand side of that estimate follows at once from the fact that sgn(a;) < if{x) 
whenever x G K. \ Si. Therefore, it remains to bound the first term on the left 
hand side of (|42|l . This is where we pick up the extra factor of (1 + A). The proof 
is essentially identical to what was done to establish (|T7|l above. Using the two 
conditions (|17|I ~ H18(I . it suffices to multiply through the following estimate by the 
quantity (1 + A)^: 

(49) II xii> \\h < II klxiV' IIl^ • II xi^.V + II \x\xi^ \\h ■ 

Here the functions xBi a-nd xBi are cutoffs which are = f on the set Bi and which 
vanish outside of 2Bi. The bound H49|) follows from evaluation of the integral: 

= / d.4x XBr ■ (^)'] dx , 



and using Cauchy-Schwartz as well as the bound |a;xBj^| ^ x'^XtSi for a suitable 
cutoff xbi ■ We have now finished our proof of the weak local energy decay estimate 



Remark 3.1. We note here that it is possible to prove (|30|l without rescaling the 
potential Q into the condition (|31|) . This can be accomplished by using the weight 
function: 

r 1 
Jo (i + e|y|)'' 

in place of H32|l above. This yields a small factor in front of |V'l(0) when it appears 
in the term on the right hand side of line (|44|l above, so one can proceed directly 
to the estimate (|^ to control things. We leave the details to the interested reader. 
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Step 3: Strong Uniform Decay of Local Energy. We are now ready to prove 
the main Morawetz estimate H19(l . With the assumptions (|13|l -l|18 |l in hand, as well 
as the weak local energy decay estimate (|3U|) . this becomes an essentially standard 
calculation. We will contract the energy-momentum tensor ((UJl with the conformal 
Killing vector-field: 

Ko ^ [t^ + x^)dt + 2txd^ = ^u^L + ^u'^L . 

The deformation tensor of this is computed to be: 

Therefore, we may form the momentum density ^^°^Pa = TapK^ and from line 
H12|l we compute the divergence: 

(50) 9" (^«)P„ = -txd^iQ) ■ - 2tQ ■ {-^f + KoW ■ H . 

By simply integrating this last line over various time slabs of the form ^ s ^ t 
and using the Cauchy-Schwartz inequality we arrive at the bound: 

(51) sup '^»)Fo(s) / / |i?(s)| • |^o(^)(s)| dxds + (^«)Po(0) 

O^s^t Jo JR 

+ [ [ [sxd^iQ) ■ (V')^ + 2sQ ■ (V^)^] dxds . 

JO JR 

Using now the identity: 

we see that (|^ in conjunction with the energy estimate implies the bound: 

(52) sup E{s) < I I \\{^ + \u\ + \n\)H{s)\\L^d.)- E^{s) ds+ m) 

O^s^t JO JR 

+ / / [sxd^iQ) ■ {ijf + 2sQ ■ [ipf] dxds. 

JO JR 

The last thing we need to do here is to bound the last term on the right hand side 
of the previous expression. We will show the bound: 

(53) ^ 

[ ( [sxd^[Q)-{i;f + 2sQ-{iljf] dxds < sn]i N^^E{s) + N{l + \^)E{Q) 

Jo JR O^s^t 

+ N f / II {l + s)il + X)H{s) \\md.)-\\ (l + A)(|Vt,,^|-f g5.|V;|) (s) \\md.) ds . 

Jo JR 

where the implicit constant is independent of the large parameter N. Notice that 
this bound substituted into immediately implies ((T^ for sufficiently large A^. 

To prove (|53|l we will chop the left hand side up into three pieces. The first is 
the "bad" set Bi. This is where most of the positivity of fL.H.S.'l H53() can be found. 
The second set is where x ^ Bi and |a;| ^ ^t. Here we use the strongly repulsive 
condition (|15|l . Finally, in the exterior of the influence of the potential when i < |a;| 
we can simply integrate things using the homogeneity bound l(TB|) . The details of 
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this procedure are as follows: 

On the set Bi we use the repulsive condition 114|) and the first term on the left 
hand side of to bound: 

[sxd^iQ] ■ (ip)^ + 2sQ ■ i^pf] dxds , 

/■* r 



{l + s)Q-{ii)Y dxds , 

'0 JBi 

< (R.R.S.)^ . 

Next, we work in the set M \ Z5i but cutoff according to how large t is. Here we 
make use of the condition (|15|l : 

1 ()rf. 

Xii^—J [sxd^iQ) ■ {il:f + 2sQ ■ {-^f] dxds , 
lOx 



JK\8i 1 + S 



< - 



(1 + s)Xi(t— )sgn(a;)a,(Q)-(V')' dxds 

JV\Bi i + S 



< (R.H.S.)IS3 • 

Finally, in the exterior where t < |a;| we use the condition (|15|l and the follow- 
ing bound which holds for parameters N such that \B2\ ^ N (where the implicit 
constant of course depends on \B2\)'- 

i-t r / IOt \ 

1 - XI (^) WdAQ) ■ W + 2sQ • (V)'] dxds , 







< 



20|52| r- f /■* 

/ 2sQ-{^f dxds ^ / / Q-{^f dxds / / Q-{tjjfdxds 
Jr Jo Js. Jn 



< N sup E{s) + N-^ sup E{s) , 



< 



(R.H.S.)® 



This completes our demonstration of (|53|1 . and hence our proof the main estimate 



4. Scattering for Small Amplitude Non-Linear Scalar Fields 

We will be brief here and leave many of the details to the reader. The main 
result of this section is the following: 

Theorem 4.1 (Scattering for Scalar Fields). Consider the Cauchy problem: 
(54) Ugcj, = A|0|f0 , m = f , dtm = 9 , 

for compactly supported functions (/, g) . Define the regularity space: 

UWn^^ = ^((l-Asp^)i 

where E_ is the Morawetz type energy from line (O . Then if 2 < p, there exists a 
universal set of positive constants £ and C depending only on p such that if: 

wmwni ^ £ , 
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then a unique solution to the problem H54|l exists for all values of the variable t and 
it obeys the bound: 

II 0(0 ^ C£ . 
In particular, one has the following uniform point-wise bounds: 

(55) \<j)\ < £ r-^ -mmil , \t-\r*\\'^ } . 

A previous result of this type was recently obtained by Dafermos and Rodnianski 
in the case of spherical symmetry and powers 3 < p (see 0)- 

To prove Theorem (j4.1(l we need four ingredients. The first is the Morawetz 
estimate (0). The second is a Poincare type estimate which will allow us to control 
the norm of our function in terms of the energy E_. The third is a paraproduct 
bound which allows us to concentrate all of our angular derivatives on a single term 
of the non-linearity A|(/)|^'(/). And the final is a global Sobolev inequality which will 
give us the bound fK5|l in terms of our energy space Hq. We now state the last 
three of these in order: 



Lemma 4.2 (Poincare type estimate for the weights u and u). Let ijj be a function 
of the variable r* , and define the weights u and u as on lines ^ . Then the following 
estimate holds: 
(56) 

j^W'dr* < [^\L4>f + uHLijf + + + - ^) (^)2 j dr* . 

Lemma 4.3 (Paraproduct bounds). On the sphere §^ the following estimates holds: 

(57) ||(l-A,,,)t(|Fm|U.(s.) < ||F||^^(g,).||(l-A,p,)tF|U.(s.) , 
for all test functions F and integers ^ k ^ p + 1. 

Lemma 4.4 (A global Sobolev inequality). Let <f> be a function of the variables 
(r*,a;). Then one has the following global bounds: 

(58) M < min{ 1 , r^(l-^)-3|i-|r*||-\ |t-|r*|pM-|kllH^ ■ 
We now give short proofs of these three Lemmas: 

Proof of estimate H56|l . The proof will follow from cutting the function %j} into three 
pieces. We write: 

V' = Xr-<-l'0 + X-I<r-<1'0 + Xr*>l'(/' , 
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where the x form a smooth partition of unity such that the dr*Xr'<±i are supported 
on the interval [—2, 2]. For the left hand portion we compute that: 



= / a,. [{t + r*)xr'<-i{^f] dr* , 

Xr'<-iW^ dr* + f {t + T*)x'r-.<^i{^f dr* +2 I {t + r*)xr'<-iM-i^ dr* 



Collecting terms and applying the Cauchy-Schwartz inequality we arrive at the 
bound: 

Xr'<~iW dr < / [i/j) dr 

+ ( f X.*<-i W dr*) ' ■( f Xr'<-i{dr'^f dr 

This easily proves (|56|l for the left hand portion of things because here one has the 
bound: 

(59) Xr'<-l{dr'^f < u^{L^f +u\L'^f . 

The proof of (|56|) for the right hand term Xr*>i4' follows from an identical argu- 
ment. We leave this to the reader. □ 



Proof of estimate H57|l . Because k is chosen to be an integer, it is easy to reduce 
this estimate to the case of Euclidean space. First of all, note that one has the 
following bound: 

ii(i-A,p,)*(|Fmiu.(s.) < iin?,(iFmiu.(s.) , 

where a is a multiindexing of the rotation generators {flij}. Via a partition of 
unity, the desired bound easily reduces to proving that for fc e N one has: 

WIF^FWh^ < II F 11^^- I! F 11^. 

Estimates of this type are well known and easy to prove (see e.g. ^J). □ 



Proof of estimate (|58|l . We first make a preliminary reduction. Because we are 
including two angular (momentum) derivatives in the norm 7Y^, via the Sobolev 
embedding on the sphere it suffices to prove the following global Sobolev estimate 
for functions ip of the variable r*: 



< min{ 1 , r2 (1 - ?^) I \t - \r* 
r ' 



(60) IV 
where H is the Hilbert Space 



\t-\r*\\~-^ } 



\n 



2 



(1 + w')(Li/^)2 + (1 + u^)iL^f + il + u^ + u^) 



(1 - 2M) 



dr* 
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Using now the Poincare estimate (|5?^ and bounds of the form l|5(^ we see that we 

have: 

(61) 



(1 - 2M) 

(i + (t-|r*|f)(a;v'f + W + {i + ij + u'y- — ^{^f 



dr* < II VII 



The first two terms on the left hand side of this last expression are enough to get 
the basic weight in the estimate (|60ll . Specifically, one has the bound: 

(l + (t~|r*|n^(V)2 < / [(l + (t-\r*\f)id:^r + W'] dr* . 

The proof of this kind of estimate is completely standard and reduces to the usual 
Sobolev bound after decomposing things via cutoff's on intervals of dyadic sizes 
according to the value of the weight (1 + (< — |7'*|)^)2. See |T| for details on this 
procedure. 

It remains for us to prove the estimate 1)60(1 for the more refined weight. It is 
clear from the form of the energy H61|) that this bound holds when r* G [—1,1]. 
Therefore we only need to establish it for the cases r* < — 1 and r* > 1. We will 
do this separately by directly evaluating some weighted integrals similar to what 
we have done many times now. For values Tq < — 1 we compute that (note that the 
weight r is essentially a constant here so we can safely disregard it): 







dr 



2Ar 1 



r 



dr* + (V.)2(0) 
r 



By collecting terms and using the Cauchy-Schwartz inequality this gives us the 
bound: 

{l-—)Ht + r*ofW'{r*o) < (l-^)^i^ (V)^(O) 
r r 

2 A'/ 1 

+ \\i^-—)Ht + r*)^j\\m[0,r^,])-\\{t + r*)dr'^\\m[0.r;]) ■ 

the crucial thing to notice here is that when the derivative falls on the weight 
(1 — '^r-){t + f*)^ it creates a positive term which can be collected with the left 
hand side. This last line gives us the desired bound because all the type norms 
are covered by the energy 7i, and so is the bound for -0 at the origin. 

To wrap things up here we need to prove the refined bound in ((60(1 for the region 
where 1 ^ r*. Notice that it suffices to do this for r* ^ t because otherwise the 
simpler weight on the right hand side of ((60(1 is favorable. Also, in this case we may 
safely ignore the weight (1 — ^) because it is essentially a non-zero constant. The 
desired bound will drop out from computing the following integral for fixed points 
Q^rl^t: 

r-\t-r*^f{^f{rl) = T dr-. [r-Ht ~ r*f{i,f] dr* + r'^ f {^,^{0) . 

Jo 
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Computing the integral on the right hand side of this last expression, throwing 
away only the term which results when the derivative falls on the weight (t — r*)^, 
collecting the terms of like sign to the left hand side, and applying the Cauchy- 
Schwartz inequality we arrive at the bound: 

r-i(t-r*f(^)2(r*) + \\ r'\t - r*)^P < r'U^ {^^{0) 

+ II r-^{t - r*)^ llL2([o,r-]) • II - r*)dr-'^p ||l2([o,^.]) . 

From this one easily derives the bound: 

r-\t-r*n^ fir*) < r'^ (V.)^(0) + || - r*R.^ |U.([o,..]) , 

which is itself bounded by the energy Ti.. This completes our proof of the estimate 
and hence our demonstration of the global Sobolev bound (jSHl- O 



Proof of Theorem \4-.l\ We now use the previous three lemmas and the main decay 
estimate (0) to prove the global regularity result. This will follow by bootstrapping 
the usual local existence theorem. We will not state or prove this local result here 
because even in this context it is an elementary application of Picard iteration and 
energy estimates (in fact, one can use the estimates developed here to set up a global 
Picard iteration because the precise structure of the non-linearity does not need to 
be preserved). Now, our theorem will follow if we can show that for each fixed time 
t up to which we have existence the weak bound supq^^^^ || 4>{s) W-hC^ ^ 2CiS implies 
the stronger bound supg^^^j || (f){s) W^-ij^ ^ C£. From the estimate 0, we see that 
we can provide this as long as we can show the bound: 

2M, 



(62) / ||(l + |ii| + |u|)(l )r.(Vl-A,p^)^(|(/.|P(/.) {s)\\L^(ar'd^) ds 

Jo ^ 

< sup II 0(s) ll^+i . 

To compute the integral on the left hand side, we use the paraproduct estimate 
(|57|l which gives us the bound: 



\\{y/l^A,p^)H\cf>\Pcl>) (s)IU.(d^) < U{s)\\l^-\\{y/l-AspH)H{s)\\LHd^) . 

Using now the definition of the energy Hq and the Poincare estimate H56|) to control 
the norm of the zero harmonic of (p, we see that we can bound: 

(L.H.S.)in3 < sup \\^{s)\\ni- f ||(l + N + |u|)(l-— )5|0|f (s)||ic»^ ds . 

The claim (I62II will now follow once we can show the fixed time bound: 

(63) \\{l + \u\ + \u\){l-—)^^\P{s)U^^, < (l + .)-i-^(P-2)||^(s)||^, . 

This last bound follows easily from the global Sobolev estimate ()58f) . To see this, 
it is convenient to split things into the three regions: 

T^i = K < -^s} , 7^2 = {|r*| Is} TZs - {^s ^ r*} . 
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In the region TZi the bound on the H63|l is completely trivial. Here the non-linear 
interaction breaks down entirely. All one has to do is to use the bound: 

(1-—) < e^"' , < , 

r 

and the fact that the are uniformly bounded (which is the best we can do!). 

In the transition region 7^2 j we use the fine bound contained on the right hand 
side of lf35|l . This gives us that: 

Taking on the weight (1 + |u| + |u|) which is ~ (1 + s) in this region, and using the 
fact that < (1 + s)-5(p-2) -^^ ^^^^^^ ^^^^ ^ 

Finally, in the Minkowski like region TZ^ we have from H58(l the uniform decay 
estimate \(t>\ < {l+t + r*)~^. This easily implies (I63|l . This completes our proof of 
the bootstrapping estimate (|62|l and hence our demonstration of Theorem 14. II □ 
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